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Abstract 

We present first principles LDA+U calculations of electronic structure and magnetic state for 
LaCo03 and H0C0O3. Low Spin to Intermediate Spin state transition was found in our calculations 
using experimental crystallographic data for both materials with a much higher transition tem- 
perature for H0C0O3, which agrees well with the experimental estimations. Low Spin state t^e^ 
(non-magnetic) to Intermediate Spin state t^e* (magnetic) transition of Co 3+ ions happens due 
to the competition between crystal field t2 9 -e 9 splitting and effective exchange interaction between 
3d spin-orbitals. We show that the difference in crystal structure parameters for H0C0O3 and 
LaCoOs due to the smaller ionic radius of Ho ion comparing with La ion results in stronger crystal 
field splitting for H0C0O3 (0.09 eV ~ 1000 K larger than for LaCoOs) and hence tip the balance 
between the Low Spin and Intermediate Spin states to the non- magnetic solution in H0C0O3. 

PACS numbers: 78.70.Dm, 71.25.Tn 
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I. INTRODUCTION 



In lanthanum cobaltite LaCoOs there is a temperature induced transition from a non- 

fl 

magnetic to a magnetic state, which has attracted considerable interest last years |l( . In the 
ground state of LaCo03 there is no magnetic moment on Co ions. At low temperature the 
magnetic susceptibility increases exponentially with temperature, exhibiting a maximum 
near 100 K. While initially this maximum was ascribed to a transition from a Low Spin 
non-magnetic ground state (tSL, S=0) to a High Spin state (t^ejj, S=2), later a new scenario 
involving Intermediate Spin state (t^e*, S=l) has been proposed. Using the results of 
LDA+U, Korotin et al. f2j explained stabilization of Intermediate Spin state over High Spin 
state due to the large hybridization between Co-3d (e g ) and 0-2p orbitals. 
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FIG. 1: Scheme of energy levels for Low Spin state of Co 3+ ions in LaCoOs. 



The physical reason for this transition is a competition between crystal-field energy A c / 
(t2 9 -e g energy splitting) and intra-atomic (Hund) exchange energy A ex (Fig. |TJ). In the 
ground state A c / is only slightly larger then A ex so that the energy of the excited mag- 
netic state of Co ion is relatively small and with increasing of temperature its population 
increases resulting in an increase of magnetic susceptibility. When La ions are substituted 
in perovskite cobaltites by other rare earth elements with smaller ionic radii, magnetic prop- 
erties show significant changes. For NdCoOs material 59 Co Knight shift measurements 3( 
showed that Co ions remain to be in the Low Spin state up to 580 K, where the grad- 
ual metal-insulator transition is observed. High-temperature diffraction study of lanthanide 
cobaltite perovskites LnCoOs (Ln=Nd, Gd, Dy and Ho) jj] led its authors to conclusion, that 
"at room temperature all cobalt ions are in the Low Spin state regardless of the Ln atomic 
number" and only as temperature increases to 1000 K "the possible electronic phase transi- 



tion can be suggested". Magnetic and electric measurements for Hoi_ x Ca x Co03 system |5( 
showed that "there is no apparent transition between Low Spin and High Spin of Co 3+ ions 
in 300-900 K temperature range" . All these data show that Low Spin state of Co 3+ ions in 
lanthanide cobaltite perovskites with smaller rare earth ions becomes more stable compar- 
ing with LaCoOs and transition temperature increases dramatically. In the present work 
we have investigated relative stability of Low Spin and Intermediate Spin states of Co 3+ 
ions in LaCoC>3 and H0C0O3 as a function of experimental crystal structure parameters 
measured at different temperatures. We have found that, indeed, while for LaCoC>3 transi- 
tion temperature value (defined as the temperature, where calculated total energy of Low 
Spin solution becomes higher than the energy of Intermediate Spin solution) is 140 K, for 
H0C0O3 this value is 1070 K. The physical reason for this is decreasing of Co-0 bond length 
due to the chemical pressure happening with substitution of La ions by smaller Ho ions and 
hence increasing of the crystal field t2 9 -e g energy splitting value thus stabilizing Low Spin 

state. Such conclusion is supported by the results of the recent diffraction experiments under 

n 

pressure for LaCoOs p, which were interpreted as "pressure-induced Intermediate-to-Low 
Spin-state transition" . 



II. CRYSTAL STRUCTURE 



LaCoOs crystallizes in the rhombohedrally distorted cubic perovskite structure 




According to the latest crystallographic data by Radaelli et al. |8[ the space group of 
this compound is R3c for all measured temperatures and has two formula units in the 
unit cell. Main structural motif of this compound is a nearly perfect CoC>6 octahedron. 
The rhombohedral distortion of the parent cubic perovskite structure can be described by 
deformation along the body diagonal so that the angle of the Co-O-Co bond has changed 
form 180° to ~163° (Fig. |2D- 

H0C0O3 has orthorhombically distorted perovskite crystal structure (prototype GdFeC>3, 
space group Pbnm) [4] (Fig. EJ), which contains four formula units in the unit cell. Like for 
LaCo03 the main crystallographic motif form CoOg octahedra, which are rotated in afe-plane 
and tilted in respect to c-axis. Due to smaller ionic radius of Ho ion comparing with La ion, 
distortion from the cubic perovskite structure is stronger for H0C0O3 than for LaCoOs and 
the angles of the Co-O-Co bond have values 149° (in afe-plane) and 152° (along c-axis). 
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FIG. 2: The rhombohedral crystal structure of LaCoO"3 (left) and the orthorhombic crystal struc- 
ture of H0C0O3 (right). Co: large spheres; O: small spheres. 



In addition to the strong bending of the Co-O-Co bond the average length of the Co-0 
bond is decreased with substitution of La by Ho (1.934 A and 1.921 A at 300 K 4, ^|). Such 
decrease can be understood as a result of chemical pressure, induced by the smaller ionic 
radius of Ho ion comparing with La ion. 

The crystal structure parameters as a function of temperature were measured for LaCo03 
in jsl and for H0C0O3 in j^J. It was demonstrated in J^J that LDA+U calculations can 
reproduce a Low Spin-Intermediate Spin transition for LaCo0 3 using temperature dependent 
crystal structure as an input. We used the same procedure for H0C0O3 in order to find if a 
small decrease of Co-0 bond length due to the chemical pressure can results in stabilization 
of the Low Spin state. 

Five 3d-orbitals are separated into triply degenerate t2 9 and doubly degenerate e g subsets 
only in cubic symmetry lattice. The rhombohedral crystal structure of LaCo03 and the 
orthorhombic crystal structure of H0C0O3 have lower symmetry than cubic as a whole but 
C0O6 octahedra are only slightly distorted in those compounds. In local coordinate system 
centered on Co ion with the axes directed to the oxygen ions t2 9 -e 5 orbitals are still well 
defined. In the following we have used such defined orbitals for analisys of the calculated 
electronic structure. 
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III. LDA BAND STRUCTURE 



The standard L(S)DA (Local (Spin) Density Approximation (t| 10]) based methods can 
not describe magnetic state for LaCo03 without taking into account Coulomb interaction 
between 3d electrons as it was demonstrated in j^j]. However analysis of the LDA results 
is instructive for understanding of the basic electronic structure of lanthanide cobaltite 
perovskites. 



We used LMTO method ll| to perform calculations of electronic structure for LaCo0 3 
and H0C0O3. The results are presented in Fig. El The radii of the muffin-tin spheres 
in LaCo0 3 for 5 K structure were RL a =3.76 a.u., Rc =2.49 a.u., and Ro=2.0 a.u. With 
temperature increasing Rc and Ro were increased proportionally to Co-0 distance for a 
given temperature. In H0C0O3 for lowest temperature (300 K) crystal structure Rc and Ro 
were taken as the corresponding values for LaCoOs scaled in proportion of Co-0 distances 
in LaCoOs and H0C0O3. Then Rc and Ro were increased proportionally to the change of 
the Co-0 distance in H0C0O3 for a given temperature. Rl & was adjusted to fill completely 
the volume of the unit cell. In the orbital basis set the following states were included: 
Co(4s,4p,3<i), 0(3s,2p), La,H°(6s,6p,5<i,4/). Partially filled 4/ states of Ho were treated as 
a pseudo-core states. 

FigureH3(left panel) shows in details LDA calculated electronic band structure of LaCo03. 
On the top of the figure total density of states (DOS) is shown. There are three distinguish- 
able sets of bands: completely filled O-p bands, partially filled Co-d bands and empty La-/ 
bands. The bands in the energy range from -7 eV to -1.5 eV originate mainly from 0-2p 
states but have a significant admixture from Co-3d states. Between eV and 2 eV there are 
empty La-4/ bands. The two groups of bands presented in the third panel, which extend 
from -1.5 eV to eV and from eV to 2.0 eV, are the Co-3<i states. Also a substantial 
0-2p contributions are apparent in this energy range. Co-3rf states, due to the octahedral 
symmetry as we discussed above, are formed by completely occupied t2 g states and empty 
e g states with a "pseudo-gap" between them. The partial densities of states for t2 9 and 
e g states are presented in the fourth and fifth panels of Fig. |SJ correspondingly. We have 
calculated centers of gravity and standard deviations of partial t2 9 and e g DOS (shown in 
Fig. EJ) in order to determine the value of crystal field splitting between t2 g and e g states. 
It was found that t2 9 states are 1.93 eV lower than e g . Also t2 9 band is much more narrow 
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FIG. 3: DOS of LaCo0 3 (left) and H0C0O3 (right) calculated within the LDA. Top panel: total 
DOS; second from top panel - partial 0-2p DOS; Last three panels: partial Co-3<i, Co-3d (tzg) and 
Co-3d (e g ) DOS. For partial Co-3d (t2g) and Co-3d (e g ) DOS centers of gravities C and standard 
deviation a of corresponding Co 3d orbitals are shown. The Fermi level is zero energy. 

than e g . 

LDA band structure obtained for H0C0O3 is shown in Fig. El (right panel). In general it 
is very similar to LaCoOs bands (except for the absence of 4/ band because Ho-4/ states 
were treated as a pseudo-core). The important difference from LaCoOa is relative positions 
of t2 9 and e g bands of Co-3d shell. As one can see from last two panels of Fig. |3] (right side), 
t2 9 states are 2.04 eV lower than e g comparing with 1.93 eV for LaCoOs and also e g band 
is more narrow (standard deviation o is 1.45 comparing with 1.59 for LaCoOs). The higher 
position of e g states and the smaller width of the corresponding band result in opening of a 
small energy gap ~0.07 eV. 

Those results can be well understood from the difference in crystal structure of two 
compounds. The values of Co-O-Co bond angles deviate from 180° much more in H0C0O3 
than in LaCoOs. That leads to a weaker d-p-d hybridization and hence to a more narrow 
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Co-3d (e g ) band. 

In adition to Co-O-Co bond bending, substitution of La ion on a smaller Ho ion lead to 
the compression of Co-0 bond lengths. Decreased value of Co-0 bond length in H0C0O3 
comparing with LaCoOa results in the increased value of hybridization strength between 
0-2p states and Co-3d states and hence to the increasing of t2 S -e 9 energy splitting of Co-3<i 
states. 

For the problem of Low Spin to Intermediate Spin state transition the most important 
is increasing of the t2 5 -e 5 energy splitting value. As magnetic transition in cobaltites is 
determined by the competition between crystal field and exchange energies, this result can 
dramatically increase excitation energy to the magnetic state and hence increase transition 
temperature for H0C0O3 comparing with LaCoC>3. 

One can estimate the change of excitation to magnetic state energy as equal to the 
change in the value of t 29 -e 9 energy splitting: 2.04-1.93=0.09 eV~1000K. This value must 
be connected to the transition temperature difference. We will show below, that this crude 
estimate agrees surprisingly well with the results of our LDA+U calculations. 



IV. LDA+U RESULTS 



In the present work we used a genera 



rotationally invariant formulation of LDA 



U ar 



proach instead of the old version used in ^ . The main idea of the LDA+U approach 3, IllJ 
is to add to the LSDA functional the term E u corresponding to the mean-field approximation 
of the Coulomb interaction in multiband Hubbard model (Hartree-Fock approximation) 

E LDA+U [p°(r), K}] = E LSDA [p°(v)} + E u [{rf}\ - E dc [{rf}\, (1) 

p CT (r) is the charge density for spin-cr electrons, E LSDA [p a {v)] is the standard LSDA (Local 
Spin-Density Approximation) functional and Coulomb interaction term E u is a functional 
of the occupation matrix n cr : 

n mm' = / I m G'inlrn,inlm'(.E)dE . (2) 

J 

Here G a , . , ,(E) = (inlma I (E — H)^ 1 I inlm a) are the elements of the Green function 

inlm,inlm v ' ' 1 v ' ' ' 

matrix in the basis of c/-orbitals | inlma) (i denotes the cite, n the main quantum number, 
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I- orbital quantum number, m- magnetic number and a- spin index) 

E u [{n*}] = 1 1 E {m} J(m,m" \ V ee \ m',m'")n" mm ,n^ m ,„+ 

(3) 

((m,m" | V ee | ml, ml") - (m,m" \ V ee \ m'" ,m!))n a mm ,n a m „ mlll }, 

where V ee are the screened Coulomb interactions among the d electrons. Finally, the last 
term in Eq. corrects for double counting and is given by 

E dc [{n a }} = l -UN{N - 1) - l -J[N\N^ - 1) + - 1)], (4) 

with N a = Tr{rtL m ,) and N = + N L . U and J are screened Coulomb and exchange 
parameters jl^ll5|. 

In addition to the usual LSDA potential, an effective single-particle potential to be used 
in the effective single-particle Hamiltonian has the form: 

H = H LSD A + ^2 I ^ mcr )^mm'(^ m '°" I) (5) 
mm' 

Km' = J2m»m»'{( m ' m " I V ee I ™' ',m">m» m "' + 

((m, m" | V ee | m',m'") — (m,m" \ V ee \ m'" , m'))n^ n ,, m w}— (6) 

U{N-\) + J{N°-\). 
The matrix elements of Coulomb interaction can be expressed in terms of complex spher- 
ical harmonics and effective Slater integrals F k as 

21 

(m, m" | V ee I m', m") = ^ a k (m, rri, m", m"')F k , (7) 

and 

A k 

a k (m,m',m",m'") = ^ (Im \ Y kq \ lm')(lm" \ Y k \ \ lm"'). 

q= — k 

For d electrons one needs F°, F 2 and F A and these can be linked to the Coulomb and 
Stoner parameters U and J obtained from the LSDA-supercell procedures via U = F° and 
J = (F 2 + F^y 14, while the ratio F 2 /F 4 is to a good accuracy a constant ~ 0.625 for the 
3d elements (l6j, |l8| (for / electrons the corresponding expression is J = (286F 2 + 195F 4 + 
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250F 6 )/6435). The Coulomb parameter U is calculated as a second derivative of the total 
energy (or the first derivative of the corresponding eigenvalue) in respect to the occupancy 
of localized orbitals of the central atom in a supercell with fixed occupancies on all other 
atoms 1^ | . 
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FIG. 4: DOS of LaCo0 3 (left) and H0C0O3 (right) calculated within the LDA+U approach for 
Low Spin state t^e^ of Co ion. Top panel: total DOS; second from top panel: partial 0-2p DOS; 
Last three panels: partial Co-3d, Co-3d fag) an d Co-3ci (e 9 ) DOS correspondingly. The Fermi 
level is zero energy. 



Coulomb interaction parameters [7=7.80 eV and J=0.99 eV used in our LDA+U cal- 
culations were computed in J^j by the constrained LDA approach j^]. The main effect of 
LDA+U potential correction (JHJ) is the energy splitting between occupied and empty states 
in such a way that the former are pushed down and the latter up comparing with LDA. In 
the result magnetic state solution, which does not exist in LSDA, becomes stable in LDA+U 
and by comparing the relative energy of two solutions, non-magnetic and magnetic, one can 
find the ground state of the system. 

We have investigated two solutions for both compounds for all crystal structure param- 
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eters corresponding to the temperatures 5-300 K for LaCoOs and 300-1080 K for H0C0O3. 
The Low Spin state solutions (Fig. 0J) are very similar to LDA results (Fig. EJ). The main 
effect of LDA+U correction is higher position of the empty Co-3d (e 5 ) states and the op- 
posite effect on the occupied Co-3d (t2 9 ) states and hence opening of a sizable energy gap 
(0.7 eV for LaCo0 3 and 1.3 eV for H0C0O3). 




Energy, eV Energy, eV 



FIG. 5: DOS of LaCoO-3 (left) and H0C0O3 (right) calculated within the LDA+U approach for 
Intermediate Spin state tJ^e* of Co ion. Top panel: total DOS; second from top panel: partial 
0-2p DOS; Last three panels: partial Co-3<i, Co-3<i (t2 9 ) and Co-3d (e 9 ) DOS correspondingly. 
The Fermi level is zero energy. 

On the other hand the Intermediate Spin state solutions (Fig. |SJ) are very different from 
LDA. At first a peak above Fermi energy appears for minority spin Co-3g? (t2 S ) states cor- 
responding to a hole in t 2g shell and at second the majority spin Co-3<i (e g ) band becomes 
half-filled (t^e* configuration). The magnetic moment values on Co ions are equal to 2.2 hb 
for LaCoOs and 2.17 \xb for H0C0O3. The Intermediate State solution gave metallic state 
in contrast to the semiconductor properties of LaCo03. In |2| this contradiction was ex- 
plained by the possibility of orbital ordering of the partially filled e 5 -orbitals of Co 3+ ions in 
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FIG. 6: Comparison of total energy per Co ion of Intermediate Spin state t^e* (full line, circles) 
and Low Spin state t^e^ (dashed line, diamonds) solutions for LaCo03 (left) and H0C0O3 (right) 
calculated with the LDA+U approach as a functions of temperature. The temperature of transition 
(calculated as the temperature where two lines cross) is ~ 140 K for LaCoOa and ~ 1070 K for 
H0C0O3. 



Intermediate State and opening an energy gap in the orbitally ordered state. Recently the 
evidence for orbital order in LaCoC>3 was found in diffraction experiments 19^ . 

The most interesting are results for total energy values of Low Spin and Intermediate Spin 
solutions. In Fig. El those energies are plotted as a function of the temperature (temperature 
is taken into account in our calculations only via experimental crystal structure parameters). 
For LaCoC>3 at 5 K the total energy difference between Low Spin and Intermediate Spin 
states is only 37 meV and (by interpolation) becomes zero at 140 K. For room temperature 
(300 K) Low Spin state is higher in energy than Intermediate Spin on 45 meV. 

For H0C0O3 the situation is very different. At room temperature there is a very large 
total energy difference between Low Spin and Intermediate Spin states (200 meV) and only 
at 1080 K Intermediate Spin state becomes lower in energy than Low Spin. The interpolated 
value of T, where energies of Low Spin and Intermediate Spin are equal is 1070 K. 

One can identify the temperatures, where lines for Low Spin and Intermediate Spin states 
cross in Fig. |H1 as magnetic transition temperatures. The value for LaCoC>3 equal to 140 K 
agrees very well with the 100 K maximum in experimental susceptibility curve [jj. The 
calculated magnetic transition temperature for H0C0O3 equal to 1070 K also agrees with 
the results of |4J, where authors conclude that "at room temperature all cobalt ions are in 
the Low Spin state" and only as temperature increases to 1000 K "the possible electronic 
phase transition can be suggested" . 
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V. CONCLUSION 



Comparative study of electronic structure and magnetic properties of LaCo03 and 
H0C0O3 was done using LDA and LDA+U approaches. The chemical pressure induced 
by substitution of La ions by smaller Ho ions leads to the bending of Co-O-Co bond angles 
and compression of Co-0 bond length that in its turn increases the value of crystal field 
splitting in H0C0O3. Total energy calculations for Low Spin and Intermediate Spin states 
demonstrated that this increase results in stabilization of non-magnetic solution in H0C0O3 
and hence increasing of transition temperature from 140 K (LaCoOs) to 1070 K (H0C0O3). 
Calculated transition temperatures agree well with the experimental estimations. 
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